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Certain intriguing consequences of the discreteness of time on the time evolution of dynamical 
systems are discussed. In the discrete-time classical mechanics proposed here, there is an arrow of 
time that follows from the fact that the replacement of the time derivative by the backward difference 
operator alone can preserve the non-negativity of the phase space density. It is seen that, even for free 
particles, all the degrees of freedom are correlated in principle. The forward evolution of functions 
of phase space variables by a finite number of time steps, in this discrete-time mechanics, depends 
^ , on the entire continuous-time history in the interval [0, oo]. In this sense, discrete time evolution 

^ ■ is nonlocal in time from a continuous-time point of view. A corresponding quantum mechanical 

' treatment is possible via the density matrix approach. The interference between non-degenerate 

quantum mechanical states decays exponentially. This decoherence is present, in principle, for all 
systems; however, it is of practical importance only in macroscopic systems, or in processes involving 
' large energy changes. 

^ ■ 

INTRODUCTION 

Time is an enigma, and many philosophers and scientists have tried to ponder over its true meaning. However, 
most physicists and mathematicians consider the concepts of space and time to be intuitively obvious, and view 
space-time as an inert and infinitely divisible continuum in which 'events' unfold. The advent of the general theory 

ly-j 1 of relativity led to the insight that space-time is dynamical and that the gravitational field should be identified with 
certain properties of the space-time continuum. Subsequent efforts at synthesizing gravity and quantum mechanics, 
either with a background space as in string theory or without a background space as in loop quantum gravity 
strongly suggest that space-time at the most fundamental level has a granular nature. In loop quantum gravity, 
the spatial Riemannian geometry is discrete |4j with the volume of space quantized in units of lp, where the Planck 

Ck length l p = (hG/c 3 ) 1 / 2 » 1.6 x 10 35 m. Further, the time evolution takes place in discrete time steps |5j of the order 

+^ ; of the Planck time t p = {HG/c 5 ) 1 / 2 « 5.4 x 1CT 44 s. 

There have been some attempts in the literature 0, & 0> to study the time evolution of Hamiltonian systems 
^ with time as a discrete parameter even before these exciting developments took place. In particular, Katayama and 
Yukawa ar gue d that, just like matter, space-time also should have an indivisible element (an 'elementary domain'). 
Yamamoto 10] also made an attempt to realize the elementary domain through quantum field theory in discrete time. 
In this paper, we propose a version of discrete-time mechanics as a precursor to the more complicated (and conceptually 
satisfying) discrete space-time mechanics. The idea is to explore the genuine consequences of discrete-time evolution, 
and not the development of numerical methods to approximate continuous-time evolution with greater precision. In 
the equation of motion approach to classical mechanics (CM), Lee ^30] found that violation of time-translational 
invariance (that accompanies discretization of time) leads to non-conservation of energy. He therefore developed 
a theory with time as a discrete dynamical variable. Subsequently, Jaroszkiewicz and coworkers |l3| succeeded in 
developing an equation of motion approach to CM with time as a discrete parameter by invoking the discrete-time 
action principle of Cadzow|14]. 

We take a fresh look at the problem and adopt the phase space density approach. We start with the Liouville 
equation for the phase space density in classical mechanics, and propose its discrete-time analogue. The structure 
of this discrete-time Liouville equation ensures that all the constants of the motion in continuous-time mechanics 
are also constants of the motion in discrete-time mechanics. The basic premises that underlie the derivation of this 
equation and the consequences of the discreteness of time (the emergence of stochasticity, correlation and non-locality) 
are discussed in Section 2. The mathematical details are relegated to Appendices A and B. We then quantize the 
discrete-time Liouville equation to obtain the time evolution of the Wigner distribution function (equivalently, the 
density matrix). It is seen that the discretization of time leads to energy super-selection and decoherence. These 
points are discussed in Section 3. Even though the present formalism is motivated by Hamiltonian mechanics, we have 
attempted to study the consequences of discrete time on more general forms of dynamical systems. In particular, the 
effect of the discreteness of time on the sensitive dependence on initial conditions of (chaotic) nonlinear dynamical 
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systems is discussed briefly in Section 4. Finally, some conclusions are drawn in Section 5. Before we proceed, we 
would like to emphasize that the fundamental unit of time r that appears in our formalism need not be equal to the 
Planck time. We also remark that the usual continuous-time mechanics is recovered in the limit r — > 0. 



CLASSICAL MECHANICS IN DISCRETE TIME 



We consider a dynamical system with I degrees of freedom described by the Hamiltonian H{x , p), where Xi and pi 
are the coordinates and canonical momenta, respectively. The time evolution of the phase space density p c t{x , p, t) 
is given by the Liouville equation 

d 

—p ct (x, p, t) = {H , p ct } = Lp ct (x, p, t) (1) 
where {H , p c t} is the Poisson bracket of H with p c t(x , p, t), and L is the Liouville operator 

'dH d dH d 

v dxi dpi dpi dxi 



J^fdH V OH n 



»=i 

Throughout this paper, the suffixes 'ct' and 'dt' as in _F ct and Fdt shall denote values of the phase space function F 
in continuous-time mechanics and discrete-time mechanics, respectively. The formal solution of Eq. Q is given by 

Pct (x,p,t) = e tL p ct (x , p , 0) (2) 

where p c t(x , p, 0) is the initial phase space density. Throughout this paper, we consider the deterministic initial 
condition 

p ct (x , p, 0) = 5{x - x (0)) 5{p- p(0))- (3) 
The value F ct {t) of the phase space function f(x,p) at time t is given by 

F c t(t) = (f(x,p)) ct = J dx J dpf(x,p)p ct (x, p, t). 

Our aim is to "deduce" the discrete-time analogue of Eq. JQ). Since there is no unique definition of a discrete time 
derivative, we do this by first stipulating the conditions to be satisfied by the equation that governs the time evolution 
of the phase space density pdt(x , p,n) at discrete time t = nr, where r is the fundamental unit of time. Specifying 
the data at one instant of time is sufficient to obtain the solution of Eq. We demand that specifying pdt(x , p ', j) 
at some discrete time t = jr should be sufficient to determine pdt(x, p, j') at any other discrete time t' = j'r. This 
implies that the equation to be discovered should essentially be a first-order difference equation in time. We thus 
obtain the most general form of the discrete-time Liouville equation for forward evolution in the form 

pdt(x,p, n + 1) - p dt (x , p, n) _ „ 

= L [ap dt (x , p, n) + ppdt[x , p, n + 1) , (4) 

T 

where < a , (3 < 1 and a +/3 — 1. Special cases of Eq. J3J are the forward and backward difference schemes obtained 
from Eq. 10} by setting a = 1 and a = 0, respectively. Restrictions on the possible values of a can be obtained by 
examining the time evolution operator T„ defined through the equation 

Pdt(x, P, n )= ( ]^_ Pdt(x,p,0)=T n pdt(x,p,0). (5) 

We assume that the Hamiltonian and the manifold in which the dynamics takes place are such that L is sufficiently 
well defined. L is then a skew-Hermitian operator, and hence its spectrum is pure imaginary. Using the relation 
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/(fi) G(L)g(fl) d?L = J G(-L)f(Q) g(Q)dn 



where 57 = (x , p) , we can show that every quantity that is conserved in the continuous time context is also conserved 
under the time evolution scheme proposed here. 
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It is clear that it is only when a = \ that T n is unitary and the evolution given by Eq. (JSJ is consistent with time 
reversal invariance. T„ is an unbounded operator (a bounded function would evolve to an unbounded one under the 
action of T n ) in the limit n — > oo for a > 5. Since this is undesirable, a must be less than 1 for n > 0. Similarly, T n 
is unbounded in the limit n — > —00 for a < i. Hence a must be greater than ^ for n < 0. Thus, an arrow of time 
emerges naturally if a ^ h. 

A more stringent constraint on the value of a can be obtained by demanding that the phase space density pdt(x , P , ri) 
be non-negative. We examine this problem in two different ways, and show that zero is the only permissible value of 
o. 



Method 1: 



Let us make a formal transformation 

(xi, ....,xi;pi, ....,pi) -> (6, ■■-,£,i;m, ,vi), 

such that L = —d/d^i . This transformation can be worked out explicitly for systems with quadratic Hamiltonians. 
While global realization of such a transformation for systems with non-quadratic Hamiltonians appears to be doubtful, 
the analysis is instructive, and we proceed with it. For the deterministic initial condition Eq. ©, the corresponding 
initial condition in the transformed variables is given by 

1 

Pdt(Cv,0) =Y[6(Zi - 6(0)) 6( Vl - 77,(0)). 
i=l 

The time evolution of pdt (£, Vi n ) is given by Eq. (J3J) to be 

Ait = 3n(£i)%i - »?i(o)) Hsfe - e 4 (o)) S( Vt - m(0)), 

i=1 

where 

As shown in Appendix A, we find that zero is the only allowed value of a which renders Pdt(£j n ) non- negative and 
T n bounded. For this case, a simple expression can be derived for g n (£i), namely, 



fc (6> - I ^ (^)'- 1 e XP (-^), taft >6(0) 
0, for < Ci(0). 



(6) 



We thus see that the backward difference scheme is the only acceptable generalization of the time derivative to the 
discrete domain. 



Method 2: 



Here we obtain a relation between p c t{x , p, t) and pdt(x , p, n) by resorting to the generating function technique. 
As shown in Appendix B, we arrive at the same conclusion: namely, that zero is the only permissible value of a that 
renders pdt(x , p, n) non-negative. This analysis leads to the following intriguing relation for n > 0: 

1 f°° 

Pdt Or, p, n) = — - I du e~ u u n ~ p ct (x , p, tu). (7) 



(n- 1)1 J 

The average value Fdt(n) of any function f(x,p) at time t — nr is given by 



1 f°° 

F dt (n) = (f(x,p)) dtin = —- ^r-J due-V-'F^ru), (8) 
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where F ct (Tu) is the phase space average in continuous-time mechanics at time t= ru. Such quantities are nonlocal 
in time from the continuous-time point of view! The discreteness of space-time thus leads to non-locality in time in 
this sense. 

Both Eqs. © and @ imply that an element of stochasticity appears in the discrete-time classical mechanics 
described in this paper. Equation suggests that the stochasticity is essentially in one variable, £1 , and that its 
conditional probability density is a gamma distribution for all Hamiltonian systems. Equation J7J) implies a random 
walk in time. Thus £1 can be identified with a kind of "internal time" . 



Examples 



We now illustrate the peculiarities of discrete-time evolution with the help of two examples: (i) a collection of free 
particles, and (ii) a collection of harmonic oscillators. 

(i) Free particles: 

For a collection of free particles described by the Hamiltonian H = Pi /(2wij) with deterministic initial conditions, 
we get the following results for the averages of Xi , Pi , xf and pf at discrete time nr: 



£>i(0)nr 



in. 



> (Pi)dt.n =Pi(0), 



<^>dt.n = + 



^(O)nrV , npl{Q)T 2 



<fi>dt,„=A 2 (°)- 



Since (pf ), n = pf (0), iJ is conserved. However, 



( x i)dt,n ~ ^)dt,„ = " 



pH0)t 2 



= Dit 



where Di = p?(0)r/m?. Hence the motion of the particles is diffusive. A similar calculation shows that 



(xiXj) 



•J/dt.n \ Xl /dt,rt \^3/dt,n 



mirrij 

which implies that even the motion of raon-interacting particles is correlated in the foregoing sense. 
(ii) Harmonic oscillators: 

Next, consider a collection of harmonic oscillators described by the Hamiltonian H — Y^i=i \ (p? + x f)> where we have 
taken all the masses and frequencies to be identical and set m j = 1 , u>i = 1 , for simplicity. The analysis of this case 
is straightforward. The final expressions for the first and second moments of X4 and pi at time nr are 



b »/dt.' 



(1 + T 2 )™/ 2 



sin (mj> + 0j(O)), 



n(0) 



(P«)dt,n= {1 + T 2 )n/ 2 cos N + 



<* 2 ) = iM 

V^/dt.ri O 



COS 



(n^ + 20i(O)) 



(1 + 4t 2 )™/ 2 



<Pi>dt,» 



r?(0) 



s (n^ + 261,(0)) 
(1 + 4t 2 )"/ 2 
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_ n (0)^(0) 

where = arctan r, (f>' = arctan(2r), rf = xf + p 2 , and ^ = arctan (xi/pi). We see that (xf) dtn - (ii)j tn and 
(Pi)dtn ~ (P*)dt n are nonzero, signalling stochasticity in and p,:. However, the Hamiltonian iJ is conserved. We 
also see that (xiXj) dt - (xi) dt (^j) dt n is nonzero, implying that all the degrees of freedom are now correlated. 



cos (0i(O)-0j(O)) 



cos (n0' + 0i(O)+0j(O))" 
(l + 4r 2 )™/ 2 



QUANTUM MECHANICS IN DISCRETE TIME 

Quantum mechanics (QM) is a remarkably successful theory with no known physical phenomena that contradict 
it. Yet, it has no universally accepted interpretation and hence the adage, "It is a theory that works for all practical 
purposes" . The founding fathers of QM have insisted that the results of measurements have to be expressed in 
classical terms. If the world is quantum mechanical, then classical mechanics (CM) should be contained within QM as 
a limiting case. Now, a crucial ingredient of QM is the principle of superposition that follows from the linearity of the 
Hilbert space. Even though manifestations of quantum mechanical interference abound in the microscopic domain, 
the application of such a principle to the macro- world seems to lead to predictions that are counter-intuitive vis-a-vis 
our day-to-day experience. The Schrodinger cat paradox, invented by Schrodinger himself, clearly brings these issues 
into focus, and it suggests that, from among the multitude of superpositions allowed by the Schro ding er equation, 
only a few robust states are allowed for macroscopic systems. The so-called decoherence program 0, 0, Hsl Hsj 
shows how this may come about entirely within a quantum mechanical description, by invoking the unavoidable 
interaction of any given system with the external world. Entanglement between the system and the environment can 
cause super-selection, i.e., the selection of a preferred set of states that are robust in spite of their immersion into 
the environment. Another consequence of the entanglement is environment-induced decoherence, which refers to the 
suppression of interference between the preferred states chosen by the super-selection rule. In what follows, we show 
how energy super-selection and decoherence may arise from the discretization of time. 

In continuous time, the density matrix p c t(t) satisfies the well known evolution equation 

9 Petit) = ^r[H,p ct (t)], 



dt 



ifi 



where H is the Hamiltonian operator and [A , B] is the commutator of the operators A and B. The usual quantization 
prescriptions whereby p c t — > p c t and the Poisson bracket {A, B} — ► [A, B]/(ih), together with the procedure for going 
from continuous to discrete time (described in this paper), leads to the following evolution equation for the discrete- 
time density matrix pdt(n)'- 



p dt (n + 1) - pdtjn) _ J_ 
t ih 



H,p dt (n + 1) =Lp dt (n + l) 



for n > 0. This implies that 



Pdt{n) = (l - tL^ (3dt(0). 



Let | a) denote the eigenfunction of H with the eigenvalue e a , and Pdt{0) the initial density matrix: 



H\a)=e a \a), Pdt(O) = / ",fflq,fl 



We then get 



Pdt{n) = a a j) 



1 + ir 



e /3 



ft 



*}</3|. 



It is clear that the diagonal elements of pdt(^) are time- invariant. The off-diagonal elements decay exponentially 
if the basis states are non-degenerate. They are, however, time invariant if the states \a) and \j3) are degenerate. 
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We thus find that, in principle, the interference involving states with different energies decays exponentially with a 
characteristic time Td given by 

T d = (2T)/\og[l + (AET/h) 2 ], 

where AE is the difference in energy between the states. If r is taken to be the Planck time (5.4 x 10~ 44 s), we find 
the decay time to be greater than 10 10 years if AE > 7meV. Thus, a microscopic system prepared in a mixed state 
by superposing states separated in energy by a few meV, decoheres only in principle, and continues to be coherent in 
practice. However, for a macroscopic system with about I0 20 particles, Td ~ 10~ 23 s for states with a 7meV change 
in the energy per particle. 

The formalism described in this paper can be extended to quantum mechanical distribution functions (QDF) as well. 
In view of the non-uniqueness associated with the classical «-> quantum correspondence (i.e., the way to construct 
quantum mechanical operators corresponding to classical phase space functions), there are many ways of defining 
QDFs 21]. We restrict ourselves to the Wigner distribution function, which follows from the Weyl correspondence 
rule. We can obtain a quantum mechanical description in discrete time by Wigner-Moyal quantization of the discrete- 
time Liouville equation. This is achieved by simply replacing the Poisson bracket by the Moyal bracket. The time 
evolution of the discrete-time Wigner distribution function Wdt(% > P,n) is given by the backward difference equation 



Wdt(x, p,n + I) - W d t(x, p,n) 



H | - sin 

h 



Wdt(x, p,n+ 1), 



where V is the operator defined through the relation 

AVB = {A, B}. 

It may be noted that Eq. (JJJ holds good in the context of quantum mechanics as well, with the re- 
placements p ct (x,p,t) — > p c t (t) and Pdt(£j P, n) — > Pdt(n) for the density matrix, and p ct {x,p,t) — > 
W c t(x , p,t) and pdt(x , p,ri) — > Wdt(% , p,n) for the Wigner distribution function. 



Connection with the Schrodinger equation 

It is interesting to ask if there exists a discrete-time Schrodinger equation whose solution is consistent with the time 
evolution of the density matrix described in this paper. If such an equation exists, then the eigenstate \a) should 
evolve to \a) n = f(n,a) \a), where f(n,a) is an unknown function that should satisfy the relation 



f(n,a)t(n,(3) 



1 + ir 



e /3 



(9) 



When a = /?, Eq. ^ reduces to \f(n, a)\ 2 — 1, which implies that fn, a) is unimodular: thus f(n, a) = exp [i 0(n, a)], 
where 0(n, a) is a real number. When this expression for f(n,a) is substituted in Eq. @, we get 



0(n, a) — ®(n, 0) — n arctan 



i~(e a - ep) 



n 

Z 7T lo S 



in which the left-hand side is real, whereas the right-hand side is complex unless r = 0. Thus the functional equation 
© does not have a solution, and hence we do not have a discrete-time Schrodinger equation consistent with the 
equation for the density matrix. This is logically consistent with the fact that the density matrix description predicts 
decoherence, which never can occur in the Schrodinger formalism for an isolated system. 



NONLINEAR DYNAMICS IN DISCRETE TIME 



We have so far concentrated on the classical and quantum mechanics of Hamiltonian systems. However, we do 
encounter more general dynamical systems while modeling numerous phenomena in physics, chemistry and biology. It 
would be of interest, therefore, to extend the prescription (for going from continuous time to discrete time evolution) 
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outlined in this paper to dynamical systems described by the set of (in general, nonlinear) ordinary differential 
equations 



~dt 



/(*) = fj({xi}), j = 1, n. 



We have the equivalent phase space formulation 

d " d 

-Q t P<A{{xi},t) = Lp ct ({xi},t) = - ^2-^r [f]{{xi\)p c t{{xi\,t)] . 

j=i o 

We note that, presented in this form, there is no obvious distinction between linear and nonlinear equations of motion 
(the differences between them would be reflected in the spectral properties of the Liouville operator). As a matter of 
fact, finite-dimensional nonlinear differential equations of motion can be reformulated as linear differential equations 
in infinite dimensions by adopting the Carleman embedding procedure|20j or similar ones. However, from a practical 
point of view, there are significant differences between linear and nonlinear problems. Of particular interest is the 
tendency of nonlinear systems to show chaotic behaviour — bounded and aperiodic evolution which shows sensitive 
dependence on initial conditions. In view of the formula (see Eq. JSJ) that relates forward evolution in discrete time to 
the entire continuous-time history in the semi-infinite interval [0,oo], it is intuitively clear that sensitive dependence 
on initial conditions in discrete-time mechanics should be different from that of continuous-time mechanics. Some 
sort of "rounding off" of the instability should occur in the discrete time context, at least in some systems that show 
bounded evolution, as the following example would suggest. 

Let the continuous- time solution for one of the dynamical variables, x, with initial value a, be given by the relation 

x c t(a, t) — cos (be ct ), with x ct (a, 0) = cos b = a, (10) 

where c > 0. The motion described by Eq. (|1(J|) is bounded for all finite t; however, it shows sensitive dependence 
on initial conditions. The ratio of the distance d c t(t) between two trajectories emanating from two infinitesimally 
close-by points (a + A) and a is given by 



d ct (t) 

to the first order in A. It is easily seen that 

dct(t) 

so that the Lyapunov exponent 



x ct (a + A,t) - a; ct (a, t) 




dx ct (a,t) 


A 




da 



|sin(&e ct )| e ct , 



lim - log (d ct (*)) = c 

is positive, signaling chaos. Using Eq. (JSJ, the corresponding expression for the distance in the discrete-time problem 
is 



ddt{n) 



d_ 

da 



1 



1 



y/T=a*(n-i)\ 



/ duu n - 1 e~ u e CTU sin(6e cr ") 
Jo 



Integrating once by parts, we get 



ddt(n) 



1 1 



1 



y/T~a7bcT (n-1)! 
Using the fact that |/ °° du cos (be CTU )f(u)\ < J °° du \f(u)\, we get 

2 1 



du cos (be CTU ) [(n - l)u"~ 



2 - u"- 1 ! e~ u 



ddt(n) < 



bcTy^T^' 
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which is bounded for all times. Thus the Lyapunov exponent for the discrete time evolution vanishes for this example. 

In order to get a feel for the differences between discrete and continuous time mechanics, let us consider a few 
more examples. If the equations of motion are such that x ct (t) <~ t a (with a > —1), then the corresponding discrete 
time evolution is given by x dt (n) ~ T a T(n + a)/(n — 1)! which goes as (iit)" for n » a. That is, a system 
that shows power law evolution in continuous time shows similar behaviour in discrete-time mechanics. If x ct (t) = 
ae bt (b > 0) in continuous time, then the corresponding evolution in discrete time is given by Xdt(n-) = ae CTn with 
c = — (1/t) log(l — for) > b. This implies that for systems that show unstable behaviour in continuous time, the 
instability is enhanced in the discrete-time context. 



CONCLUSIONS 



To conclude, we have shown that the backward difference scheme preserves the non-negative character of the phase 
space density of a classical Hamiltonian system. The time evolution operator is a bounded operator for all times. We 
have shown that discretization of time leads to stochasticity. Irrespective of the number of degrees of freedom, one 
function of the phase space variables becomes stochastic. This variable has a unique probability density which turns 
out to be the gamma density. The discrete-time evolution by a finite amount in the forward direction depends on the 
entire forward-time history of the continuous-time evolution. In this sense, the discrete-time evolution is nonlocal in 
time. The formula that relates Pdt(x , p,n) to p ct (x , p, t) is suggestive of a random walk in an internal time. The 
motion of even free particles becomes correlated motion, in our formalism. 

The same formalism is amenable to a quantum mechanical treatment via density matrices (equivalently, via the 
Wigncr distribution functions). The elements of the density matrix that connect degenerate states are time-invariant, 
whereas the ones that connect non-degenerate states decay exponentially. Thus discretization of time leads to energy 
super-selection and decoherence in quantum mechanics. It is also interesting to note that an arrow of time emerges in 
the present framework. Finally our approach does not seem to permit a consistent description of quantum mechanics 
via the discrete-time Schrodinger equation. 

Some of the consequences of the discretization of time that we have pointed out in the foregoing are intriguing, and 
warrant further investigation to elucidate their meaning and implications. The complexities that may arise when the 
Liouville operator is not sufficiently ' 'good should also be explored. 

To conclude, classical stochasticity, quantum decoherence and an arrow of time emerge automatically in the present 
version of discrete-time mechanics. All these features disappear and we recover the usual classical and quantum 
mechanics in the limit r — > 0. 
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CALCULATION OF g n (^) 
We now proceed to calculate g n (£,i) given by the relation 

*«->-(£^)" «(*-<■<»»■ 

The Fourier transform g n (k) of g n (Ci) is given by 



Rewriting (1 + iak)/(l — i/3k) as (3 — ij3rk) 1 — a] and using the binomial theorem, we get 

«*)-(3)"tG") ( -^(™) 1 ^ ro - 



Taking the inverse Fourier transform of Eq. , we get 

9n(ti) = (-f) n s(ti-m) + (J)"e (") (12) 

where hj(^i) is given by the relation 

h(n = ljrw&- ^(o))'" 1 ex P h 15 ^ 1 ]. fOT & > MO) 

\o, fore,<ei(o). 

In Eq. tlH , the first term is singular, while the second is the sum of a finite number of regular functions. The first 
term is negative for odd n. It then follows that a has to be zero for the probability density g n {£i) to be positive for 
all n. In that case g n (£i) = h n (£i). 

TIME EVOLUTION OF p dt 

For forward evolution, i.e., n > 0, define the generating function 

G(f, p, z) = ^z n p dt (x, p , n) = z™ ( J />dt(x, P, 0), 

n=0 n=0 ^ P T J 

to get 

G(x , p, z) = ( ) p dt (^, P, 0). 

\ 1 — z — (p + azjrL J 

This expression can be further simplified to 

G(x, p, z) = ( — — ) p dt (x, p, 0) + ( — - — ) ( ^— — J p dt (x, P, 0). 

\p + az J \p + az J \ 1 — z — \p + azjrL J 



Now use the integral representation 

1 



1 — z — (/3 + olz)tL 



r-OO 

I du exp {— u [1 — z — {/3 + az)rL]} 
Jq 

y * (i + aTL ) m / d u 



u ra e -u e liTL 



10 



to get 



G(x,p,z) = f^^p dt (x,p,0)+(-^-)f2^(l + arLr 
\ p + az I \ a + az ) ' m! 

v ' ' m— 

x / du u m e~ u p ct (x , p , (3tu) , (13) 



where p c t(x , p, t) is the continuous-time phase space density given by Eq. (J2J). Collecting together the coefficients of 
z n yields 



- x ( ay _ _ l^V «Y J (l + a^' 
p dt (a;,p,n) = p dt (a; , p , 0) + - ^ J 

/>00 

x / duu e~ u p ct (x , p, (3tu). (14) 
Jo 

The first term on the right-hand side in Eq. (|14fl is singular, and negative for odd n. The second is a sum of 
finite-order derivatives of a regular function. Therefore, in order for pdt{x , p, n) to be non- negative, a must vanish. 
Thus, both the methods referred to in the main text lead to the same answer. However, the second method yields 
the representation 

1 f°° 

pdt{x, p, n) = -. — r / due~ u w" -1 p ct {x , p, ru) 

[n - 1)! Jo 



for p dt (x , P, n) in terms of p ct (x , p, t). 



